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Abstract 

This paper investigates disentanglement as a result of evolution according to a class of master 
equations which include dissipation and interparticle interactions. Generalizing an earlier result of 
Diosi, the time taken for complete disentanglement is calculated (i.e. for disentanglement from any 
other system). The dynamics of two harmonically coupled oscillators is solved in order to study 
the competing effects of environmental noise and interparticle coupling on disentanglement. An 
argument based on separability conditions for gaussian states is used to arrive at a set of conditions 
on the couplings sufficient for all initial states to disentangle for good after a finite time. 
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1. INTRODUCTION 



Entanglement - the fact that to all density matrices can be written as 



P 1&2 = $>/^ 2 (1) 

i 

with pi a probability distribution - has become one of the most notable differences between 
classical and quantum mechanics. What at first sight seems a rather technical difference, 
lies underneath all 'no local hidden-variables' theorems which help set quantum mechanics 
so firmly apart from classical mechanics (Refjlj, RefjsJ). More recently, entanglement has 



come to be regarded as a resource on which all essentially quantum technologies are based 
(quantum computation, cryptography, e.t.c, see e.g. Ref|3j). An understanding of how en- 
tanglement comes to be destroyed by interaction with an environment is therefore important 
not only to explaining how the classical world emerges from an underlying quantum one, but 
also to understanding and controlling the effects of noise, which make quantum technologies 
such a practical challenge. 

In Refpj, a theorem was proved categorizing completely disentangling maps. We call a 
completely positive trace preserving map 



M ■ B{H) -> B{H) (2) 



completely disentangling if 



M ® ljv : B{H ® C N ) -> B(H <g> C*) (3) 

has M®ljv(p 1&2 ) disentangled for all density matrices p 1&2 G B(H<g)C N ) and for all positive 
integers N. Ruskai Ref[(| showed that M is completely disentangling if, and only if 

M(p) = ^2tr{fi k p)m k (4) 

k 

with {/ife} a POVM and m k density matrices. 

Diosi has taken advantage of this to show that certain markovian open-system evolutions 
result in a map of this form after a certain finite time Ref P| . This finite time environment- 
induced effect at first seems quite strange when compared with the more familiar asymptotic 
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environment-induced decoherence effects of similar time-scale. But it should be noted that 
separable density matrices form a convex, finite-volume region within the space of all den- 
sity matrices. If the equilibrium state is separable, we expect the evolution to 'suck' states 
towards it along smooth trajectories which, at a certain time, pass into the separable region 
around the equilibrium state. This picture makes clearer how disentanglement occurs in 
finite time, whereas decoherence, where states are 'sucked' smoothly towards a stable sub- 
manifold of measure zero (the states diagonal in the relevant basis), can only be achieved 
asymptotically. 

In Refj^| Diosi and Kiefer showed that the wigner function evolution 

dW _ p dW D d 2 W 
dt m dq pp dp 2 

results in a positive P function after a finite time. Diosi used essentially the same technology 
to show that this evolution is also completely disentangling after a finite time. The positive P 
function result itself implies a disentanglement result if we consider the obvious extension of 
(JHJ to an iV-particle system: a positive P function is a classical probability distribution over 
the (separable) coherent-state projectors. The complete disentanglement result if stronger 
however, and also gives disentanglement on the symmetric iV-particle system. 



This paper is an elaboration and continuation of Refllol. and is designed to be read 
in conjunction with it. In Ref[16], the relation between decoherence and disentanglement 
was discussed and a new technique, based on separability criteria for gaussian states, used 
to compute the finite-time disentanglement of general states in a bipartite system, when 
evolved according to (j5J). The evolution of an E.P.R.-like state was solved and used to 
illustrate how the separability criteria come to be satisfied. In this paper we first focus on 
extending Diosi's argument regarding the complete disentanglement of a single system to 
more general, dissipative evolutions. In particular we will be interested in evolutions which 
are of Lindblad form, i.e. 

^ = -i[H, + { 2L kP L t " LkLlp ~ P L i L >) (6) 

n 

with non-zero dissipation. Then we will also go on to apply the technique of Ref[16| to 
a harmonically coupled pair of particles, where there is competition between the coupling, 
which can generate entanglement, and the environmental influence, which tends to destroy 
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it. We will concentrate on determining conditions which guarantee that the environment 
wins. 

The layout of this paper is as follows. Below, in Section 2, we extend Diosi's complete dis- 
entanglement result, analyzing a more general and realistic evolution. For a specific choice 
of Lindblad operator, we compute the complete disentanglement timescale for high temper- 
atures relative to the dissipation. In Section 3, we then go on to consider the open system 
evolution of bipartite systems where the particles are coupled. Complete disentanglement 
results cannot be immediately applied to this situation because the interaction will tend to 
generate new entanglement between subsystems. We will use a different method, based on 
the recent separability criteria for gaussian states, which explicitly produces a decomposition 
of the state into separable states. In section 4, we conclude and make some remarks about 
the role of entanglement destruction in explaining classical behaviour. 



2. DIOSI ANALYSIS FOR A GENERAL DISSIPATIVE MASTER EQN 

We start with the following Markovian (quantum Fokker-Planck) master equation for the 
Wigner function 



dW pdW 2 dW n d(pW) „ d 2 W 
-77— = 7: h mu) q— h 27 



dt m dq "* dp ' ' dp ^ qq dqdq ^ pp dpdp ' qp dqdp 

This describes a harmonic oscillator interacting with an environment whose effects can be 
well approximated by white noise. In practice, this is true fairly generically of commonly 
encountered environments at high temperatures. 

n 

It can be shown (see e.g. Ref|10|) that this master equation is of Lindblad form if, and 
only if, the diffusion matrix D, defined by 



d 2 W d 2 W 

+ 2D av ^Z- (7) 



D 



is positive and 



^ qq ^ q%> 
Dqp Dpp 



(8) 



det£> > 



7 

4m 2 



(9) 
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In analysing this master equation, we will make much use of the Fourier transform 

W = / ^-±-W{z)e i{ *^ (10) 



where 



(a,b) = a T b (11) 

is the standard euclidean inner product, 



(12) 




and 




and similarly without 'bars'. 

The Fourier transform W obeys the equation 



dt m dq dp dq 

with 



which we can solve by characteristics. 



(13) 



dW pdW 2 dW n dW T - - 

- - qmu 2 —- + 2jq—- - z T DzW (14) 



D= I DpP DqP I (15) 



2.1. Specialization to the free case, to = 

Diosi's argument is given in Ref[16(. Here, we explain the argument as we go, for u = 0, 
but considering 7^0. The calculations for the case w ^ are contained implicitly in the 
section below on coupled systems. We do not pursue Diosi's analysis in detail for this case 
because, as noted below, oscillating terms destroy the monotonicity of disentanglement and 
therefore make an analytical treatment much harder. 



5 



The gist of the argument is to show that, by smearing things out in phase-space, the 
propagator for the evolution quickly results in a positive P function. This ability of open- 
system dynamics to smear out the negative parts of phase-space quasi-distributions has 
been studied before with regard to treating them as genuine probability distributions in 
pseudoclassical circumstances. Here, the crucial observation is that a map which always 
results in states with positive P function is a completely disentangling map, of the form 
All this is accessible because the propagators are gaussian in form, meaning the calculations 
can be carried through explicitly, and because the properties of gaussian states are well 
understood. 

We start with the evolution written in the form 

dW p dW dW T - - . N 

— = ^— + 2jq— - - z T DzW 16 
ot m oq oq 

which is equivalent to 



dW(z t ,t) 



dt 

where the equations for the characteristics are 



zfDz t W(z t ,t) (17) 



dqt Pt „ _ 

-77 = 2 7 g t 

at m 

^ = (18) 
dt v ' 

Since this set of equations is linear and Markovian we can write 

z t = E t z (19) 

where E t is a family of matrices affecting the evolution as a 1-parameter group. We easily 
find 





(20) 



Setting 

IH= I drEjNE T (21) 
Jo 
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we then can solve as 



W(z,t) = W{E t {E^z),t) 

= W(z- t ,0)exp[-(z- t ,ntz-t)} (22) 

Taking the inverse Fourier transform, for the original Wigner function we have 



W(z,t) = J ^W(z^,0)eM-^-t^tZ-t)]e l{z ^ 

= J J rf 2 z , e- i ^''^-* > ^(/,0)e i ^' , ' z " ) exp[-(^ t ,/x^_ t )]e i<z '^ ) 

= J ' d 2 z' J ' ^IdetiE^expit^Ej^-z^^jexpi-iz^tZ^Woiz') (23) 

We will not be so interested in the evolution of the argument, nor the normalizing prefactors, 
which take care of themselves. The middle Fourier transform can be carried out to give the 
propagator as 

/ (^2 exp ^ z,r? ^ exp-(z,^ t z) 
= J -0^exp-(z + i^^r]z,/2 t (z + i^^r]z))e-^ z ^ lriz) 

27Tvdet fit 4 

~ g(2v T »tv) (24) 

The notation g(E) denoting a normalised gaussian with covariance matrix E. In short then 

W t (z) = g(2r] T pi t rj) * W (v T Efrjz)\ det E t \ (25) 

* meaning convolution. 

We calculate 2rj T [itV to be given by 



2r) T m = 2 



n ( 4 7 t+4e- 2 ^-e- 4 ^-3 ^ , p p 2~ l t+ e - 2 ~< t -1 , p , n 2 7 t+ e - 2 ^-l , p f 

2 7 f+e- 2 ^-l , D . D , 

^PP 47 2 m ^ -^ppf- 



(26) 
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We will now need the fact that the P function and Wigner function are related by 



W(z) = 5 (C 1/4 ) * P{z) (27) 

with 

I 1 / D pp rn I 



2 V 2 

Because the convolution of gaussians factors, (where it makes sense) we have 

P t {z) = g{2r] T m - C lfi ) * W^r? Ej r)z)\ det E t \ (29) 
and as for Wigner functions 

W p * W a = tr(por) > (30) 

the P function will be > if g(2r] T fi t V — C1/4) is a Wigner function. In 2 dimensions this is 
true if, and only if, (RefQ) 

det(2r7 T /i^-^i/4) > i (31) 

This positivity of the P function manifests the density matrix as a probability distribution 
over coherent state projectors, which is therefore of the form (j3J). 

For long times (i.e. once the exponentials and constants are negligible), we have 



, 7 ,/ \A>Pl67 3 m 2 + 2-D gP4 2jra + D qq t j [D pp 2m +D qq t 

2r] /i t r] - Ci/4 ~ 2 v ' t , ' v 7 I (32) 

Due to cancelations, the determinant is 

2 2 det(D)t 2 > (33) 

provided the master equation is Lindblad (by (JUJ)). Thus it is clear that we have complete 
disentanglement in finite time. Exact determination of the timescale requires numerical 
methods. 



S 



2.1.1. The limit of no dissipation: 7 — > 

The limit 7 — > allows comparison with Diosi's calculation. One finds 

3m2 m 2m (34) 

^ I 

which, allowing for differences in notation, is the same as in Refp and gives (for a master 
equation of the form (J5|L i.e. for D qq = D qp = 0) the disentanglement timescale as calculated 
there: 



2.2. A particular choice of Lindblad operator 

We consider the 'minimally invasive' modification that results in an evolution equation 
of Lindblad form, and which satisfies the fluctuation-dissipation theorem (see e.g.Ref^j]). 
This means we take the Lindblad operator to be 



L = \fAmkTx + ij\ — p (36) 

V AmkT 



resulting in a master equation with 



D pq = 

D pp = ImkT^ 

D <« = s^r (37) 

This equation has been used in the study of realistic models, such as the quantum brownian 
motion model in certain regimes (e.g. Ref[12]). Diosi's criterion 

det(2r ? r ^7 7 - Cyi) > - (38) 
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becomes 



det 



9n 4 7 f+4e- 27t -e-^-Z , on , 1 



(^-Dpp 7 2 7 2 m + 2D qq t 2 



reducing to 



> 



^Dppt 



Dppin 



1 

4 

(39) 



$ x ( (-^) 2 ( 7 t) - ) + ( 7 t)< 




;^i)2 e 2 + ( ^i )0 > o (40) 

7 7 



or, better, 



___ (4r _ i) _ (2 + l( 7 t ) 2 )r + ( 7 t )V 



1 

4' 



e 



(7*o) 4 (7^o) ; 



> 



(41) 



with 



(42) 



and 



$ = 4( 7 t )r + 4e~ 27 *° T - e" 7 * " - 3 

= 2( 7 t )r + e- 27i0T - 1 (43) 
and where we have introduced a decoherence timescale, to, given by 

t^VW < 44 > 

The decoherence timescale - the time taken for reductions of e -1 in quantities measuring the 
progress of decoherence - is significant because, for reasons discussed in Ref jl(| and as found 
in practice, disentanglement and decoherence occur on comparable timescales. Decoherence 
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is gene rally more familiar and has been better studied than entanglement breaking (however 
sec 

RefQ, RefQ), and so it is of interest to compare the two timescales by using the ratio 

T. 

If we are at 'high' temperatures, i.e. with respect to the damping (or put differently, 
if decoherence occurs much more quickly than damping, as is the case for situations of 
emergent classicality) : 

7*o = ^ « 1 (45) 

we can conjecture a solution to the above equation in the form of a power series in (7^0) ■ 
To second order we have 



$ = ( 7 to)(-3r) + ( 7 to) 2 (r 2 ) (46) 
9 = ( 7 t ) 2 (2r 2 ) (47) 

which yields the complete disentanglement time r*(as a fraction of the decoherence time) 

r. = \ - fhio) 2 (48) 

This then represents the timescale for this system to become completely disentangled, i.e. 
disentangled from all other isolated systems. This particular form shows that damping does 
not affect the overall behaviour, and in this regime shortens the timescale for its occurrence. 
The result is not directly comparable with (|35|) because of the distinct choice of D. 

If the damping is extremely strong, a better result may be obtained by using (J3l3|) directly 
to estimate the timescale. 



3. DISENTANGLEMENT OF COUPLED SYSTEMS 

In this section we consider the separation of two harmonically coupled particles, each in 
interaction with an environment. Without the coupling the particles might evolve according 
to (say) 

dW _ pi dW d 2 W _ p 2 dW d 2 W 
dt mi dqi pp dp\ m 2 dq 2 PP dpi 
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The 'reduced' Wigner function for particle 1 would then evolve according to exactly 
By the disentanglement results above, the particle-1 Wigner function would have to be 
completely disentangled after a finite time. In particular, this means that the full Wigner 
function for particles 1 and 2 separates. As observed above, the results on complete disen- 
tanglement specialise to results on iV-partite separation, unless the particles are coupled in 
some way. 

In this section therefore, we consider the case of two particles which are allowed to 'talk' 
to each other via a harmonic potential. (More general potentials being considerably more 
difficult to treat.) A potential can act to build up correlations between each party and can 
countermand some of the disentangling effects of environmental noise. 

Once we allow the particles oscillating modes, we must abandon all hope of simple analytic 
expressions estimating disentanglement times. Graphs of quantities describing the level of 
entanglement now will have 'wobbles' in them, in addition to any simple overall trends. The 
precise disentanglement times will depend in detail on the relationship between parameter 
values, which determine the size and phase of the wobbles. What is more, there is no reason 
to think that the coupling could not allow the pair of particles to resurrect their entanglement 
(in fact, see Refjl^). We therefore content ourselves with determining conditions on the 
couplings that guarantee disentanglement between the pair for good, after a finite time; 
without worrying about the precise details and timings of the separation. 

3.1. The Evolution 

We will use the 'capital' , rotated coordinate system given by 



Qi 



l 



(Qi + Q2) 



1 ~ 7! 



Q2 



1 



(?i - 92) 



2 = 



1 



(Pi + V2) 



P 



1 ~ 7! 




1 



(50) 
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or 



h -h 

The differential equations relevant to the method of characteristics are 



—rr = 2 7<2i 

dt m 

—7- = 2 7 g 2 

dt m 
dPi 

~dT 

dP 2 
dt 



z := Rz (51) 



= mQ, Qi 

= (2mu 2 + mtt 2 )Q 2 := mVL l2 Q 2 (52) 



where u is the frequency of the inter-particle coupling and Q the strength of the well. 
The 1 and 2 evolutions then decouple and take the same form: 



dt 



(53) 



^ —27 — mT l 
^ muj 2 

This decoupling will allow us to pursue most of the calculations in the 'dashed' and 'un- 
dashed' sectors in 2 dimensions, and in parallel. Where we use the same notation at both 
the 2x2 and 4x4 level, we intend it to be understood that the 4x4 matrix is block diagonal 
and made up from the dashed and undashed matrices of the same name. 
Solving this differential equation we get 

, e~ 7 *(cosat + ^ sin at) _ e ^tsmat 
Et = ma (54) 

mVl 2 e - ltS m^t e --rt(cosaf +2sincrf) ' 



where a = \/Vt 2 — 7 2 

The diffusion matrix remains the same in the rotated coordinates. We will take the 
diffusion coefficient matrix to be diagonal, having in mind eventually to use the minimal 
Lindblad values as above. We obtain 
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(55) 



where 



a = T I 1 + + + — sm 2ar + cos 2ar lT (1 " ^ " 

o = — — 1 sm Zar + cos 2ar 

\2ma 2 2a 2 J \ 7 J 

2ri 2 ^ 2 T\ + — sm2ar + oos2ar (^-(1 - -) - ^-jj 

(56) 

and 



/£>!, 7 2 x Vfi 4 \l-^ 21 ' D17 , ^ 2 ^ 

* = (t (1 + ^ ) + ^)^^ + 2^T^y( a - Qe 1 '-2 at - 7 e-^ S1 „2 Qt ) 

(D 1 7 2 D 2 m 2 fi 4 ^ 7-7e- 27 *cos2at + ae- 27 *sin2at 

B = 



a 2 ' 2a 2 J 2( 7 2 + a 2 ) 

/ Di7 D 2 mfi 2 7\ /e~ 27 * - 1 a(a - ae" 2 ^ cos2at - 7c" 27 * sin 2at)\ 
\2ma 2 2H 2 J V 27 2 7 ( 7 2 + a 2 ) J 



C = (J?l— + £l(l + £-)) ± e27 ' + D o 21 - (a -ae- 2 ^cos2at - 7e -^ sin 2at) 
V2m 2 a 2 2 v a 2 7 2 7 2a( 7 2 + a 2 ) v ' 7 

7 2 ^ £>i \ 7-7e~ 2 ^cos2at + ae- 27i sin2at 

+ ~ ~ 2^V J 2( 7 2 + a 2 ) ( ' 

With all this, we have 

W t (Z) = g Mt * W (e t Z)\ det E t \ (58) 

where 

et = r) T E?r) (59) 
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and 



(60) 



3.2. Showing separation 

Now we have solved for the evolution, how are we to address the issue of separation? 
Simple, general criteria for separability of states are notoriously difficult to come by. There 
are three notable exceptions: H = C 2 £g> C 2 ; 7i = C 2 <8> C 3 ; and, embedded within Ti. = 
L 2 {R) g) L 2 (M), the class of gaussian states. Separation criteria for gaussian states come 
in a variety of flavours (see RefS and for relations between them RefQ); we shall use the 
particularly simple one from the paper of Duan et al.. In terms of our rotated coordinate 
system, this criterion reads: a Wigner function, g(£), is separable if, and only if 



Separation criteria for gaussian states will be useful to us because the propagators for our 
evolution are gaussian. In words: if , after a certain time, we can take off part of the 
propagator and use it to smear the initial Wigner function positive, and still be left over 
with a gaussian that separates, our evolution will be a separating one. In symbols: 



Sn + S 44 > 1 



E22 + £33 > 1 



(61) 



W t {Z) 



gu t * W (e t Z) 



detE, 



■ e* t (g Mt *W )(Z) 



(9e- t M t eT_ t *e* t W )(Z) 
(9e- t M t eT t -C *9C* ^W )(Z) 



(62) 



which means we can examine to see whether 



1. the Gaussian g e _ t M t e T t -c separates 



2. we can simultaneously choose C such that gc * ^Wq > 0. 
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If we can do both these things, then W t (z) separates. It will be of the form 



5> fl y d 2 z[d*z' 2 p(z[, z' 2 )WM( Zl - z[)WP(z 2 - 4) (63) 



where 



p = g c *e* t W \detE t \ (64) 
and the (not necessarily discrete) sum over a describes the separation of the state g € _ t M t e T t -c- 

3.3. Calculations relevant to separation 

Pursuing our second point, it is straightforward to see that 

g c *e* t W = e* t (g etCe T*W ) (65) 

so that this will be > if g et c e f is a Wigner function. When is g ttCt T a Wigner function ? 
The 'capital' or rotated coordinate system will be used unless otherwise stated. If we have 
a gaussian g M with 

( m \ 

M = (66) 
\ m! J 

and the matrices m and m' give rise to 2 x 2 gaussians each of which is a Wigner function, 
then g M is Wigner. We know that this is the case if 



i 

m H — i] > 

w! + -r] > (67) 
2 

which in 2 dimensions is equivalent to 



1 

det m > - 
4 

det to' > - (68) 



as long as the traces are positive. 
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In our chosen coordinate system, the relevant matrices may be taken to be block diagonal. 
The majority of the calculations then follow along in 2 dimensions. We pick C block diagonal 
and such that 

e t CeJ = C (69) 
with Co positive enough to give a Wigner function, i.e. 

det Co = - (70) 
4 

in both the dashed and undashed sectors. Thus g et c e j is a wigner function by construction. 

Consistent with this, we must now address our first point. We want g e _ t M t e T t -c to be 
a separable Wigner function. First, we must check it is a Wigner function. This will turn 
out to be the difficult condition to satisfy, in the sense that the separation criterion, while 
messy to check, follows on trivially if g e _ t M t e T t -c does become a Wigner function. We work 
a the level of the 2x2 blocks again 

We want 



det(e. t M t e T _ t - C) 
= det E 2 _ t det{2 fi t - (r]C r] T )) 

= e 4 ^det(2/i t -(^ r ))>i (71) 



We choose 



. my/^kT 1/2 \ 

vCoV 1 = \ ' (72) 

1/2 l/2m^kT J 



and for reasons explained above, declare ourselves only interested in the asymptotics, and 
conditions guaranteeing finite-time termination of entanglement. Therefore we throw away 
any term containing e -27 *. We also proceed with D\ and D 2 chosen as in the 'minimal' 
Lindblad case above. This means 
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D l( ^ + 37 2 ^-4 7 4 ) ggn^g 

4 7 fi 2 (fi 2 - 7 2 ) + 4 7 + 1 J 
mfcT(^ 4 + 3 7 2 ^ 2 - 4 7 4 ) m^ 2 _ 2y 
2ft 2 (fi 2 - 7 2 ) 32kT [C ' 

7 fcT7 , ^,__ 27t . 



16A;T ft 2 



+ 0(e" 



r Dl i A^ 4 + 3 7 2 ^ 2 - 4 7 4 ) 

4m 2 7 fi 2 4 7 fi 2 (fi 2 - 7 2 ) 1 ; 

= ^-+ (^ + 37 2 ^-4 7 4 ) 27t 

2mft 2 32mA;TO 2 (ft 2 - 7 2 ) 1 J 1 J 



So, ignoring terms 0(e 



, ^ Tu (Q 4 + 3 7 2 fi 2 — 4 7 4 ) / (kT) 2 , kT n 2 , 



(ft 2 - 7 2 )ft 2 I ft 2 V 7 (16A;T) 2 ' V kT 



(g + 3 7 2 ft 2 - 4 7 4 ) 2 _ ft 2 _ _ JT7 2 

+ 7 + 16ft 4 (ft 2 - 7 2 ) 2 32v^Tp n 2 [ Q 2 } 

/ 7 x 2 , 7 2 2A;T 7 7 



, +- - + — ^ (74) 

v 8fcT' 2ft 2 ft 2 8fcT K ' 

If this last quantity (consisting of various dimensionless comparisons between the constants 
present) is positive definite, our Wigner condition will be satisfied on a timescale set by 7 . 
It is not transparent which regimes this will be positive in, but at high temperatures, the 
important term is 



2 ^" 4 + 3 7 2 ft 2 - 4 7 4 ) x ^ (ft4 + 3 7 2 ft 2 - 4 7 4 ) 

7 2 (ft 2 ~ 7 2 ) J ^ osc / ldec > \ 7 2( fi 2_ 7 2) 

with t dec = to and t osc = ft -1 . This will be positive if the system is under-damped, i.e. 

ft 2 > 7 2 . 

Note: if the particles are both sitting in a harmonic well, this condition will need to be 
satisfied for both ft and ft'. 

We must now check that separability will also occur. Above. (|61|). we stated the Duan et 
al. criteria for separability. Note that one could equally use r}T>ri T instead of E. Applying 
this to g t _ t M t t T ' -Ci we become interested in 
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V e_ t (M t - C )e T _ t rj T = E£ t (2fit - (rjC V T ))E- t 

Now 



(Ql,2 + P 2,l)9 E T t ^ t _ (vCov T ))E _ 



[ dZ(Q\ 2 + P^) - 1 =exp 

/ dX{E T _ t X) T M 1 , 2 {E T _ t X)g 2 ^ t 

-(r?Co?? T 



2El t (2(it - (vCor] 



{{E T _ t X) T mi , 2 {E T _ t X)) (2x2) +{{E T _;X') T m 2 , l {E T _;X')) c ) 



= tv{E T _ tmi , 2 E_ t (XX T ) (2x2) )+tr(E^ t 'm 1>2 E_ t '(X'X T ') (2x2) ) 

y 2 M -( v C 0V T ) y 2 M J-(^C ( ' jJ7 T ) 

= tr( J E T ,m 1 , 2 E_ i (2^ - frCWO)) + ti(Ef£ t 'm lt2 E'_ t (2^ - (vC' oV T ))) 
= ■ (1,2) + (1,2)' 



where 



M 1)2 



mi = 



m 2 = 



and everything is still block diagonal. 




E T _ t m x E_ t = e 27 * 

and 



(cos at — 1 sin at) 2 (cos at — 1 sin at) 

v a ' ma v a ' 

^ (cos at - 2 sin at) (^) 2 

ma v a ' \ ma / 



E T _ t m 2 E_ t = e 



2~ft 



/ mgWt N2 - m " 2 sinaf (COS at - ^ sin at) 

v a / a v a ' 

mil 2 sin at 7 ^, + \ 7 ™vA2 



sm at (cos at — 1 sin at) (cos at — 1 sin atY 

a v a / v a ' 



19 



In fact 



tx{E T _ t m li2 E_ t N) = vl 2 Nv lt2 (81) 

where are the vectors given by the first and second columns of E-t read upwards. This 
means that 



in A- ; (2 /f/ (,}CW))) = e 2 *u T ( 1 ° J (2^ - ( V C r] T )) f 1 7/< * | „ 

-7/a 1/ma I \ 1/ma 



(82) 



and 



1 1 \ „ / — mfi 2 /« 

tr(^V2^- i (2/i i -(r/C r ? T ))) = e 2 ^u T | (2^-(^ )) | // 

— mfi 2 /a —7/a y \ 1 —7/a 

(83) 

where w T = (cos at, sin at), a euclidean unit vector. This is still a mess, but this form 
allows us to see that, as long as we satisfy (f74"j) > (the Wigner condition), the separability 
condition will come to be satisfied on a timescale set most importantly by 7. 



3.4. The case of 7 = 

We now consider the case of 7 = 0. For 7 = and a diagonal diffusion matrix, we can 
proceed as above with much simplified equations. We find 



j D pp D qq m 2 Vl 2 \^ + | f D pp D qq m 2 Vl 2 \ sin fit 



2 2 J \ 2 2 J 2tt 

B = 

c ~ {2^ + ~) t+ {~~2^n 2 )^r (84) 

The requisite conditions will obviously come to be satisfied after a finite time, guaranteed 
by D > 0: the propagator spreads out in a linear way and the evolution is 



cos fit 



_ sin Qt 
mil 



E t =| mu I (85) 

mfi sin fit cos at 

which clearly doesn't stop E'^_ t /j, t E t from coming to satisfy the separability criteria. 
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4. DISCUSSION 



Given the similarity between the results in Section 2 and Section 3 it is probably worth 
summarising and discussing in what ways they are related. In Section 2, we considered the 
evolution of a system according to (J7J). We showed that after a finite time, the map effecting 
the evolution came into the form 

M{p) = ^ti{p k p)m k (86) 

k 

which meant that system had become disentangled from all other systems. What happens 
if the system of interest is itself composite? If we imagine that the p in (JSfiJl is actually a 
state on A ® B, then we could enquire about the form of the reduced state on A after such 
an evolution. Let us expand the initial state p as 

ij 

where the {vr^} are not necessarily positive. Then we have 

tT B M(p) = 22 tlA ( tr -s (w)) tr £* m fc 

k 

= ^tr A (pfTTij (ti B Pkpf)) tr B m k (88) 

kij 

The map governing the evolution of trgp will therefore not usually be of the completely 
disentangling form unless p is factored in the first place. 

An exception to this is if the systems A and B do not interact, in the sense that their 
evolution map takes the factored form M = Ma®M B - This is the case for instance, when the 
coupling to of Section 3 vanishes, as in (j4T?j) . When this happens there are closed equations 
governing the evolution of each subsystem, and they each result in a completely disentangling 
evolution. In particular then, systems A and B must disentangle. Seen in this context, the 
point of Section 3 was to investigate how much coupling between two systems spoils their 
disentanglement by a noisy environment. The answer to the question is essentially that it 
does not, at least for the conditions of interest for emergent classicality. 

Of course, Section 3 also illustrated how the technology devised to test for separability in 
gaussian states can be used to investigate dynamical disentanglement in bipartite systems 
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for quite general states. 

All this has further confirmed expectations that uncontrolled environmental interactions 
are highly destructive to entanglement in quantum systems. That this destruction occurs 
(and occurs quickly) is not surprising given our experience of decoherence in such systems, 
and since entanglement relies on the superposition principle for its existence. It is interesting 
though, that by contrast with decoherence, entanglement - even with all isolated systems 
- can be wiped out totally after a finite time. When the entangled systems are allowed to 
interact with each other, the situation is less simple as the interaction can help to rebuild 
lost entanglement. However, in certain regimes (i.e. for certain regions in the space of 
temperature, damping coefficients e.t.c), the entanglement between the systems will still be 
totally destroyed after a finite time. 

Although it may be of practical interest to know about the effects of environmental in- 
teractions on entangled systems (so as to avoid it in quantum technologies), in investigating 
these questions we have been primarily interested in emergent classical behaviour. Entan- 
glement is an inherently quantum phenomenon, and so it is encouraging to find an efficient 
mechanism by which it is rapidly destroyed in the systems of our everyday experience. But 
to what extent, and in what sense, is the absence of entanglement necessary for a system to 
admit a classical description? Entanglement between two systems A and B certainly means 
that we cannot hope to describe their correlations in a classical way. For example, if they 
remain entangled under evolution, we would not expect to find decoherent histories which 
refer to each system separately. In the situations where we normally try to explain the exis- 
tence of a classical description however, there is not normally this clean division into definite 
subsystems. Given a gas of particles in a box, what sort of entanglement properties do we 
need to destroy in order to be able to describe it classically? There is no quantitative answer 
to this as yet. But a picture hinted at here is that we would like the subsystems which we 
refer to by the coarse-grained observables of classical interest to disentangle. For instance, 
suppose we were interested in local number, energy and momentum densities: n(x), h(x) 
and p(x)(see e.g. Ref jl^l] ) . Near equilibrium, we expect to be able to describe the system 
by states which are approximate eigenvectors of n, h and p. What we mean is that there 
should be no entanglement with respect to these indices in order for us to be able to find a 
classical description. 
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